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The Scale of Inflation in the Landscape 
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We determine the frequency of regions of small-field inflation in the Wigner landscape as an 
approximation to random supergravities/type IIB flux compactifications. We show that small- 
field inflation occurs exponentially more often than large-field inflation The power of primordial 
gravitational waves from inflation is generically tied to the scale of inflation. For small-field models 
this is below observational reach. However, we find small-field inflation to be dominated by the 
highest inflationary energy scales compatible with a sub-Planckian field range. Hence, we expect a 
typical tensor-to-scalar ratio r ~ C2(10 -3 ) currently undetectable in upcoming CMB measurements. 
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I. INTRODUCTION 

Recent years have seen both the advent of precision 
cosmology giving strong indications [1-9] for an early 
phase of cosmological inflation [10-14], and theoretical 
evidence for an exponentially large landscape of meta- 
stable dc Sitter vacua [15-20] combined with the first 
models of inflation in string theory [21-23]. As the num- 
ber of inflationary model realizations and final states pro- 
vided by dS vacua with small vacuum energy is quite pos- 
sibly extremely large, a description of inflationary observ- 
ables is in need for a statistical description if one wishes 
to move beyond the lamp posts given by existing model 
constructions. 

Inflationary models are generically sensitive to the 
presence of higher-dimension operators (e.g. from ra- 
diative corrections or integrating out heavy fields), and 
this sensitivity naturally splits the model space into two 
parts [14]. In small-field models of inflation [11, 12] the ef- 
fective canonically normalized inflaton scalar field evolves 
parametrically less than a Planck distance in field space 
during the 60 efolds of cosmologically necessary inflation- 
ary expansion. Control of dimension-six corrections to 
the scalar potential is sufficient for this class. Large-field 
models [13] involve the inflaton crossing a parametrically 
super-Planckian distance A^o during the same 60 efolds. 
In such models, successful slow-roll inflation necessitates 
the suppression of corrections at any dimension which 
amounts to the presence of a protecting symmetry [14]. 
The only extant symmetry capable of protecting large- 
field inflation and which has been embedded into string 
theory so far has been a shift symmetry of an axion- 
like pseudo-scalar field. These axions arise generically 
in string compactifications [24-27] where they can yield 
large-field inflation using monodromy [28]. 

Generically, these two classes are accompanied by 
an observational discriminator. Inflation produces pri- 
mordial curvature perturbations and gravitational waves 
with nearly scale-invariant power spectra (A^ ~ H 2 /e, 
and A j- ~ H 2 , respectively) originating as quantum fluc- 
tuations stretched to super-horizon wavelengths. The 
fractional power in gravity waves (tensor modes) r = 
Aj-/ ' A\ = 16e is controlled by the first slow-roll param- 
eter e = Ckin/^H 2 < 1. Its smallness enforces a vacuum- 
energy like equation of state during inflation which is nec- 



essary to drive accelerated expansion. For a large class 
of models the slow-roll of the inflaton translates into a 
monotonically increasing evolution of e. This leads to a 
relation between A^60 and the scale of inflation H which 
implies that large-field inflation is necessary to produce 
a sizable tensor mode fraction r > 0.01 in reach techno- 
logically during the next few years [29]. [52] 

By being tied to the scale of inflation, the tensor mode 
fraction r is an inflationary observable which will at most 
have a statistical description on the landscape. Hence, we 
need to determine the distribution of inflationary vacuum 
energies for accessible regions of the landscape. A guid- 
ing motivation here is that an analysis of the distribution 
of extremely small vacuum energies close to zero on the 
landscape has already been successful in providing an an- 
thropic explanation of the smallness of the observed pos- 
itive late-time cosmological constant (c.c.) [15, 35]. The 
vacuum energy distribution very roughly factors into a 
contribution coming from a number count of inflationary 
solutions, and a cosmological factor which involves vac- 
uum transitions described by tunneling events [36] and 
the subtleties of eternal inflation. 

Recent work has analyzed the cosmological probability 
distribution factor [37]. This led to the surprising answer 
that the physics of tunneling-mediated vacuum transi- 
tions and eternal inflation largely decouple from the dis- 
tribution of vacuum energies parametrically smaller than 
the Planck density. Hence, the cosmological prior is flat 
which leaves the inflationary vacuum energy distribution 
on the landscape to be determined to leading order by 
model realization and vacuum counting. We are thus 
left with comparing the relative number frequencies of 
small-field and large-field inflation models on an acces- 
sible region of the landscape which we here choose to 
be the landscape of type IIB flux compactifications on 
warped Calabi-Yau manifolds (CYs). 

Hence, in this note we determine the number frequency 
count of small-field inflation models on the landscape of 
supersymmetric type IIB CY flux vacua. Using random 
matrix theory, we find that there are exponentially many 
more small-field inflation models in the moduli potential 
of the type IIB flux landscape than there are proper dS 
vacua. Comparing this with the restrictions on large-field 
models occurring on this landscape discussed in [37], we 
therefore statistically expect the absence of primordial 
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tensor modes r > 0.01 in upcoming CMB observations. 

II. THE WIGNER ENSEMBLE AND RANDOM 
SUPERGRAVITIES 

The F-term potential of Af = 1 supergravity 

V = e K (F A F A - 3\W\ 2 ) (1) 

is the starting point of the analysis of critical points in 
the landscape. As usual Fa = 3aW + W8aK and W 
and K are the superpotential and the Kahler potential 
respectively. Critical points are defined by the condition 

d A V\ cp = (2) 

and can be maxima, minima or saddles. To determine 
the nature of a given critical point one must analyse the 
eigenvalues of the Hessian matrix, defined in terms of 
the F-term potential as H mn = d mn V where m,n can 
be holomorphic or anti-holomorphic indices. Taking into 
account the structure of the F-term potential of Eq. (1), 
the Hessian decomposes into a sum of the form 

H = HsiJSY + H K (3) +% p ure + T~L K (i) +V. s hift- (3) 

s , ' > . ' 

W ishart+W ishart Wigner 

Each of these matrices is defined in terms of the Kahler 
potential, the superpotential and their derivatives [19, 
38] . For our purposes it suffices to review the definitions 
and some properties of the Wishart and Wigner matrices 
(for a review see [41-43]). 

A Wishart matrix [40] is a complex matrix defined as 
M = AA^ where A is a random Nf x Nf complex matrix 
drawn from some distribution with mean [i and variance 
a: Q(n,a). Its eigenvalue spectrum has support on the 
interval [0, 4iV/<r 2 [, is peaked towards the origin and is 
given by the Marcenko-Pastur law [44] . 

A Wigner matrix is a Hcrmitian matrix defined as M = 
A + A\ where A is drawn from a distribution J7(/z,cr). 
The eigenvalue spectrum of the Wigner ensemble is given 
by the Wigner semi-circle law 

which can be obtained by unconstrained integration of 
the joint probability density function (pdf) 

dP(\ 1 ,...,\ Nf ) = exp(-±J^\A n(A,-A,) 2 (5) 
\ i=i / i<j 

over all but one variable. Equation (5) gives the probabil- 
ity of generating a matrix with eigenvalues in [Aj, A; + 6 A] 
and it will be crucial for the analysis of the probability 
of inflation in the landscape of random supcrgravities we 
will present later. A rather useful physical interpretation 
of Eq. (5) was put forward by Dyson in [45] in terms of 



a one dimensional gas of charged particles moving under 
the influence of an attractive quadratic potential and a 
repulsive mutual interaction. This picture proves very 
useful in qualitatively estimating behaviour of the sys- 
tem. 

A crucial property of the eigenvalue spectrum of the 
Wigner ensemble is that for the cases of interest, in 
which the random matrices are drawn form a distri- 
bution 0(0, l/s/2Nf), it has support on the interval 
[— y/2, \^2] Mp. So unlike the Wishart ensemble, which 
has all eigenvalues positive, a typical Nf x Nf matrix in 
the Wigner ensemble will have Nf/2 tachyonic directions. 

The typical eigenvalue spectrum of random supcrgrav- 
ities, as defined by "H, was found analytically in [38] 
through the free convolution of the constituent spec- 
tra. The spectrum has support in ~ [—0.7, 7.5] Mp 
(for Minkowski vacua) and so it typically features sev- 
eral tachyonic directions, meaning that the most likely 
critical points in random supergravity are steep saddles 
rather than a local minima. 

While the eigenvalue spectrum of the full random su- 
pergravity is distinct from that of a Wigner matrix, it is 
certainly true that its tachyonic part has its origin in the 
Wigner matrix since the spectrum of the sum of Wishart 
matrices is positive definite. 

The presence of the positive semi-definite contribution 
from the Wishart matrices in the full random supergrav- 
ity leads to a substantially enhanced frequency of local 
minima compared to a Wigner matrix based estimated. 
However, as the frequency of inflationary regions relative 
to local minima is dominated by the tachyonic part of 
the spectrum originating in the Wigner matrix spectrum 
alone, this relative likelihood of inflation is still deter- 
mined to leading order by the Wigner matrix estimate in 
the full random supergravity as well. Conversely, the ab- 
solute frequency of inflationary regions will be enhanced 
in the full random supergravity proportional to the in- 
creased occurrence of local minima. 

Studies of the string landscape often involve compu- 
tation of the probability of occurrence of critical points, 
with particular emphasis on minima, suited for descrip- 
tion of the present day Universe. These spectra corre- 
spond a large the shift of the smallest eigenvalue to the 
right of its typical position and are exponentially unlikely 
[38, 39, 46]: 

Pmin ~ e- cN * + °W p„ (l). (6) 

In this letter we analyse small field inflation in the 
same light and try to determine how likely it is to find 
sufficiently flat saddle points in the landscape using the 
Wigner ensemble as our main tool. The reasons to ap- 
proximate the full Hessian by a single Wigner matrix are 
twofold: firstly it is the Wigner matrix that gives rise to 
the tachyonic directions and so by focusing on these one 
hopes to uncover the inflationary structure behind the 
full Hessian; secondly for the Wigner ensemble we are 
in possession of the joint pdf, Eq. (5), whose numerical 
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integration allows us to estimate probabilities without re- 
curring to direct counting. The joint pdf that lies behind 
the full Hessian of random supergravitics, Eq. (3), is 
unknown and so direct counting, the generation of large 
samples of matrices and the counting of the ones that 
have the spectra we are looking for, is the only probe 
available. Since we are looking for minima and flat saddle 
points, which are extremely rare events, direct counting 
is computationally expensive. 

We therefore focus our analysis on the Wigner ensem- 
ble, presenting the results in the next section. 



III. INFLATION IN THE LANDSCAPE 

We start by deriving an identity regarding the proba- 
bility for inflation in the Wigner landscape. As explained 
above, the distribution of saddle points in a random su- 
pergravity will be given by the Wigner ensemble as the 
leading approximation to the full supergravity Hessian. 
By simple manipulation of the integration limits it is pos- 
sible to prove that inflationary saddle points are exponen- 
tially more abundant than minima with masses greater 
than the inflationary mass. For our purposes, g-field in- 
flation happens in a saddle point in which q fields have 
masses in the range [—77,77] and Nf — q fields in [77, oo[, 
for suitably small 77 > 0. 

The probability for generating a Wigner matrix with 
all eigenvalues greater than —77 can be found by integra- 
tion of the joint pdf: 



Nf 00 

P(VA >-77) = II / dAidP(Ai, ■■•> Ajv>) 

II J dA 4 J] J dAjdP. 



E 

Tl = 



n\(N f - 77)! 11 



(7) 



In going from the first to the second line of (7) we 
have simply split the integration region into [— 77,00 [= 
[—77, 77[U[t7, oo[ for each A, taking care to include the cor- 
rect combinatorial factors. Using Dean and Majumdar's 
result regarding the probability of large fluctuations of 
extreme eigenvalues for the Wigner ensemble [46] 



P(VA>£) = e- 2 *«K, 
where $(£) is given by 



1 

108 



36^ 



one may write Eq. (7) as 
P(inf) 



27 log 18-2 log(-£ + y/6 + e) 



(8) 



(9) 



P(VA > 77) 



e 2AcN s 



(10) 



with Ac = $(77) — $(—77). Henceforth P(inf) denotes the 
total probability for inflation, defined as the sum over all 
possible inflationary dynamics for a given Nf, i.e. 



Nf 



(11) 



In a manifestation that it is statistically more expen- 
sive to displace the lowest eigenvalue to 77 than to —77, we 
see that Ac > and so flat saddle points, suited for in- 
flation, are exponentially more frequent in the landscape 
than minima with all masses larger than 77. 

The main aim of this work is to determine the ratio 
P(inf)/P(min), where we define P(min) = P(VA > 0). 
Once again the results of [46] allow us to push ahead. 
Noting that 



P(min) 
P(VA > 77) 

one finds 

P(inf) 



,-2(*(0)-*(77))JVj = -2AcNj 



(12) 



P(min) 



(e 2AcW| _ 1)e 2AcAT| _ e 2^'(0)N] + 0{j] 2y 



(13) 

We therefore expect inflationary saddle points to be ex- 
ponentially more abundant than local minima in the 
Wigner landscape. 

In order to confirm and extend the above results we 
estimate the relevant probabilities by Montecarlo inte- 
gration of Eq. (5), setting 77 = 0.1, in the window 
Nf 6 [2,16]. We then fit the relevant probabilities for 
each value of Nf to the exponential law of Eq. (8) as is 
expected from the- theory of large eigenvalue fluctuations 
developed in [46]. The results are presented in table I. 
We see that our method systematically overestimates the 





Analytical 


Fit 


P(\ > -r,) 


0.447 


0.429 ± 0.004 


P(min) 


0.549 


0.530 ± 0.004 


P(A > 77) 


0.665 


0.645 ± 0.004 


P(inf) 




0.403 ± 0.002 



Table I: Analytical estimates and fits to numerical data. 

probabilities of occurrence of these rare events. This is 
reflected on a shift of the fitted parameters on the level of 
a few percent. We stress that even though the error bars 
cannot account for this deviation, the fact that the nu- 
merical and analytical results show the same trend lends 
credibility to our results. 

In Fig. 1 we plot the probability for finding an infla- 
tionary saddle point in the landscape, presenting both 
the data points, the analytical estimate [47] 



P(inf) = e~ 2 ^ N f - e~ 2 ^ N f . 



(14) 



and the best fit of Table I. 
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Figure 1: Probability for inflation as a function of Nf. In 
blue (lower line) the analytical estimate of Eq. (14) and in 
red (upper line) the best fit of Table I. 

As anticipated flat saddle points, like minima, are ex- 
tremely unlikely in the Wigncr landscape as they cor- 
respond to large fluctuations of the smallest eigenvalue. 
However since it is statistically costlier to displace the 
smallest eigenvalue to than to — r\ = —0.1, flat sad- 
dle points are exponentially more abundant than local 
minima as is illustrated in Fig. 2. The ratio given by 

Pjinf) j e QA27N f fitted 

P(min) ~ \ e °- W9N f analytical ' 1 j 




l k— -1 

2 4 6 8 10 12 14 16 

Nf 

Figure 2: P(inf)/P(min): Inflationary saddle points are ex- 
ponentially more likely than local minima in the Wigner land- 
scape. In blue (lower line) the analytical estimate and in red 
(upper line) the best fit. 

We will now relate this behaviour in terms of the cutoff 
77 on the mass of the fields to the 2nd slow-roll condition 
Vv = v/V < 1- For this purpose, we note that our re- 
sults above were obtained by choosing the variance of 
the Wigner ensemble to be a = 1/ y/2Nj. This approxi- 
mates a random supergravity where the mass eigenvalues 
distribute according to the Wigner semi-circle law on a 
range [—y/2, y/2] in units of Mp. The crucial point to 
observe is that a typical supergravity landscape has both 



its typical potential energy and mass eigenvalue scale 
characterised by the gravitino mass 777,3/2 = e K ^ 2 W as 
this controls the typical size of the individual contribu- 
tions in (1): |(V)| - m 2 3/2 ~ y/ ((didjV) 2 ) . Therefore, 

the choice a = 1/ ^J2Nf with its typical mass eigen- 
value size of 0(1) describes random supcrgravities with 
777,3/2 ~ Since for such supcrgravities we then 

also have \(V)\ ~ m %/2 ~ we have 77 ~ ijy and 

a cutoff 77 < 1 in the integrations of (7) directly im- 
plies slow-roll. The study of actual string theory de- 
rived example landscapes [16, 48-50] points to scenar- 
ios where \(V)\ ~ mj /2 < M£ UT - 10" 5 . We can 
now use the Wigner semi-circle law (4) together with 
the joint pdf (5) to rescale a — > am,y 2 which will ap- 
proximate the mass eigenvalue distribution of a random 
supergravity with |(^)| ~ \J ((didjV) 2 ) ~ m \/2 an< ^ 
eigenvalue range [-^"13/2,^^3/2]- This forces us to 
rescale the integration limits in (7) to zLr/m 2 ^. As we 

now have a/ ((didjV) 2 ) ~ m^ 2 , we now get that the 2nd 

slow- roll parameter r\y = r\ m 2 / 2 /\/ (i^i^j^) 2 ) ~ 77 is 
again specified by the original cutoff 77 < 1. Therefore, 
the exponential enhancement which we found above for 
7773/2 ~ 0(1) generalises to the known string landscape 
regions which can be approximated by random super- 
gravities with TO3/2 Mgut controlling both the typical 
size of the scalar potential and the mass matrix eigen- 
value size. 

Note that this exponential enhancement is estimated 
conservatively, as the random matrix description of the 
critical points of a random supergravity by definition 
selects for either minima or saddle points. Yet, small- 
field inflationary regions do exist on almost flat inflection 
points of the scalar potential as well, with a tuning cost 
comparable to that of flat saddle point. Therefore, our 
method is conservative in that it underestimates the to- 
tal rate of small-field inflationary regions occurring in a 
given random supergravity. 

The same method that lead us to the above conclusions 
also allows us to discern what is the preferred inflation- 
ary dynamics for a given Nf. Dyson's interpretation of 
Eq. (5) in terms of a gas of charged particles gives us 
a hint of what behaviour to expect. For any particular 
value of Nf there are Nf possible types of inflationary 
dynamics: from single field to Nf field inflation. Single 
field inflation corresponds to having only one eigenvalue 
in the range [—77,77] and the remaining Nf — 1 in [77, 00 [. 
For large values of Nf this is highly unlikely since eigen- 
value repulsion in the interval [77, 00 [ would tend to push 
one or more eigenvalues into the inflationary region. On 
the other hand Nf field inflation is also very rare, since 
it corresponds to squeezing all eigenvalues in the narrow 
range [—77,77], leading to a configuration where the re- 
pulsive force would tend to push some eigenvalues out 
of this interval. Somewhere between these two limiting 
cases one can find the most likely behaviour. In Fig. 3 
we plot the ratio P(q — inf)/P(inf) as a function of Nf 
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for rj = 0.1. 
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Figure 3: Relative likelihood of q-field inflation as a function 
of the total dimensionality of the field space Nf. 

We observe that the transition from single to two field 
inflation happens at Nf = 5 with the next transitions 
from 2 to 3 and 3 to 4 field inflation happening around 
8 and 12 respectively. We note that the values for which 
the various transitions happen depend strongly on rj: the 
larger the rj the sooner the transitions will happen. A 
quantitative understanding may be developed by study- 
ing the distribution of spacings between adjacent eigen- 
values. 

Next, we recall that the minimum total number of e- 
folds of slow-roll inflation at a critical point scales with 
?y as Ntot ~ l/M [47, 51]. The question of whether we 
should select for the maximum amount of slow-roll in- 
flation (due to the maximised 3- volume growth) or not 
amounts to a choice of the measure of eternal inflation. 
Therefore the answer to the question whether we expect 
single-small-field or multi-small-field inflation to domi- 
nate the small-field regime likely depends on the choice 
of the measure. 

The presence of several fields contributing to inflation 
close to a saddle point or inflection point has the poten- 
tial of generating local non-Gaussianity which is absent 
in the single-field case. As this is tied to the relative 
importance of single-field versus multi-field, statements 
about possible non-Gaussianity emanating from a multi- 
small-ficld regime again likely depend on the choice of 
the measure. We leave this for future work. 



IV. CONCLUSIONS 

In this note we have determine the number frequency 
count of small-field inflation models on the landscape of 
supersymmetric type IIB CY flux vacua. As the effective 



4D theory of both dS vacua and small-field modular infla- 
tion models in this region of the landscape is described 
by 4D Af = 1 supergravity, we have used random ma- 
trix theory to describe the region's vacuum structure in 
terms of a random supergravity [19, 38]. Meta-stable dS 
vacua require a fully positive-definite mass matrix (Hes- 
sian). Such Hessians constitute an exponentially sup- 
pressed fluctuation of all eigenvalues to positivity in the 
context of the theory of the random Hessians from 4D 
M = 1 supergravity. Consequently, we expected small- 
field inflation models which relax the positivity for at 
least one of the eigenvalues to be favoured compared to 
full mcta-stability. Our analysis of the Wigner ensem- 
ble giving the leading order description of this effect in 
random supergravity matched this expectation. We find 
that there are exponentially many more small-field infla- 
tion models in the moduli potential of the type IIB flux 
landscape than there are proper dS vacua. The analysis 
of the frequency of large-field models and the cosmolog- 
ical probability factor in [37] led to an estimate for the 
relative likelihood of large-field inflation 

Pa 06O >m p ^ {hbl) fry>i (16) 

P&<t>60<M P Pflat saddle 

We may now plug in that [37] (3 h i.i >1 < 1 (not all 
CYs will support the topological requirements for ax- 
ion monodromy) and (h2 ) < h 1 ' 1 < 0(100), as well 
as our results here f3 f i at saddle ~ exp(+ScNj) > 1. Fi- 
nally, we note that the small-field model enhancement 
Pfiat saddle is the largest for the least tuned saddle points 
with | ^7 1 ~ O(0. 1). Upon imposing COBE normalisa- 
tion on the generation of inflationary curvature pertur- 
bations, these saddle points also have the largest energy 
scales of all small-field models, and in turn only mod- 
erately sub-Planckian field ranges. Hence, we predict 
small-field inflation to dominate abundantly, and to be 
concentrated at the largest energy scales compatible with 
a sub-Planckian field range. Consequently, we expect a 
typical tensor-to-scalar ratio r ~ 0(1O -3 ) which may be 
within reach of future CMB B-mode polarisation mea- 
surements. 
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